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Abstract 



We compute the free energy density for gauge theories, with fermions, 
at high temperature and zero chemical potentiaL Specificahy, we analyticahy 
compute the free energy through 0{g'^), which requires the evaluation of three- 
loop diagrams. This computation extends our previous result for pure gauge 
QCD. 
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I. INTRODUCTION 



The perturbative expansion of the free energy of high-temperature gauge theory has the 
form 

F ~ T'[co + c^g' + csg' + {c'Jng + c,)g'' + 0(/)] , (1.1) 

where the Cj are numerical coefficients (with some dependence on the choice of renormahza- 
tion scale) and where we have assumed the temperature high enough that fermion masses 
can be ignored. In a previous work we showed how to compute the coefficient C4 of 
g^ in pure, non-abelian gauge theory from three-loop diagrams. We shall now incorporate 
fermions into the theory and so obtain a three-loop result for QED and real QCD. This 
computation is a mostly straightforward extension of our previous work, and so we refer the 
reader to that work for motivation and pedagogy. In fact, the basic calculations we need to 
do for fermions very closely parallel those we did previously for bosons, and our object in 
this paper will simply be to point out the relevant differences, catalog results for the basic 
building blocks of 3-loop calculations, and present our final results. 

In the next section, we fix our notation and conventions for coupling constants, group 
factors, and so forth. In section 2, we outline the basic integrals that are needed in order 
to compute fermionic contributions to the free energy. In section 3, we show how to derive 
analytic results for those basic integrals, though many of the details are left for appendices. 
Finally, in section 4 we present our result for the free energy and discuss its sensitivity to 
the choice of renormalization scale. 

II. NOTATION AND CONVENTIONS 

We'll consider gauge theories given by classical Euclidean Lagrangians of the form 

Ce = ^P {d, - igA^T-) 7^^ + \ {d,Al - d^A^ + gP'^AlAl)' + (gauge fixing) , (2.1) 

where the T° are the generators of a single, simple Lie group, such as U(l) or SU(3). To 
simplify presentation, we will not derive results for an arbitrary product of simple Lie groups 

2 



such as SU(2)xU(l), but such cases could easily be handled simply by adjusting the overall 
group and coupling factors on the results we give for individual diagrams, and Ca are 
the dimension and quadratic Casimir of the adjoint representation, with Ca given by 

dp is the dimension of the total fermion representation [e.g. 18 for six-flavor QCD), and S-p 
and 5'2F are defined in terms of the generators T" for the total fermion representation as 

Sp = -^tr(T2) , = ^i^iT'f] , (2.3) 

where = T°'T'^. For SU(N) with Uf fermions in the fundamental representation, the 
standard normalization of the coupling gives 

dA = A^'-l, Ca = A^, dp = Nn,, Sp^-m, 52F = ^^nf. (2.4) 
For U(l) theory, relabel g as e and let the charges of the Uf fermions be g^e. Then 

dA = l, Ca = 0, dF = nf, SF = J2^i^ S2F = Yl^i- (2-5) 

i i 

We shall work in Feynman gauge. We also work exclusively in the Euclidean (imaginary 
time) formulation of thermal field theory. We shall conventionally refer to four-momenta 

— * 

with capital letters K and to their components with lower-case letters: K — {ko,k). All 
four-momenta are Euclidean with discrete frequencies ko = 27mT for bosons and ghosts and 
ko = 27r T for fermions. We regularize the theory by working in c? = 4— 2e dimensions 

with the modified minimal subtraction (MS) scheme, which corresponds to doing minimal 
subtraction (MS) and then changing the MS scale // to the MS scale fi by the substitution 

To denote summation over discrete loop frequencies and integration over loop three- 
momenta, wc use the short-hand notation 



PO 



-n 



[IV 
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•n;v = O 



(b) 



FIG. 1. The (a) bosonic and (b) fermionic contributions to the one- loop gluon self-energy. 



for bosonic momenta and 



p 



(2^)3-2^ 



for fermionic momenta, where 



E 

{po} 



E 



P0=27r(n+i)T 



We shall also sometimes use the notation 



(2.8) 



(2.9) 



P T{P} 



(2.10) 



rp+{p} 

We handle the resummation of hard thermal loops (which is required to make perturba- 
tion theory well-behaved beyond 0{g'^)) as we did in ref. IQ. Specifically, we must improve 
our propagators by incorporating the Debye screening mass M for Aq, which is determined 
at leading order by the self-energy diagrams of fig. 0: 



,2 rab 



M's" = nS(o) = n" (0) = g's 



This is accomplished by rewriting our Lagrangian density, in frequency space, as 



f2.111 



(2.12) 



where Sp^ is shorthand for the the Kronecker delta function Sp^^. Then we absorb the first Aq 
term into our unperturbed Lagrangian Cq and treat the second term as a perturbation. 

III. THE BASIC INTEGRALS 



The most basic one-loop integrals that appear in high-temperature field theory are of 
the form 



bn = 



fn — 



1 



(3.1) 



,p p2n ' •' " jf^pj^ p2n ' 

The bosonic form of these integrals needed for the calculation were reviewed in ref. ^ and 



are given by 



\ ( . n C'(-l) \ 

+ 0U). 



- + 21n-^ + 27. 



(3.2) 
(3.3) 



As noted in ref. 0, the /„ are then easily determined by considering fn + bn and then scaling 
the momenta {po,p) by 2 so that it becomes proportional to 6„. One finds 



fn = {2 



2n+l-d 



l)bn. 



(3.4) 



In ref. [0], we reviewed how three-loop diagrams contributing to the free energy can be 
reduced to some simple sum-integrals at 0{g'^). The most basic was 

1 



rbb 



(3.5) 



ipQR P^Q^R\P + Q + Rf 
which corresponds to the basketball diagram of scalar theory, depicted in fig. 0(a). 11^ is 



defined by 



(3.6) 
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(a) (b) (c) 

FIG. 2. The (a) bosonic, (b) fermionic, and (c) mixed scalar basketball diagrams. Solid (dotted) lines 
correspond to scalar propagators with bosonic (fermionic) frequencies. 





• • • • 



(a) (b) 

FIG. 3. The (a) bosonic and (b) fermionic setting sun diagrams. 

and we have introduced the superscript b for 11 to indicate that it is defined with a bosonic 
frequency sum. When fermions are included, the reduction of diagrams to a few simple 
integrals requires introducing some fermionic relatives of /bail- 

lLii^f^[^'iP)]\ (3.7) 
I^in^f^ll\P)Il'{P), (3.8) 

These are depicted by figs, ^b) and (c). 

Another basic integral encountered in the pure gauge theory case was the one associated 
with the scalar sunset diagram of fig. |5(a), evaluated to leading order in masses: 



When fermions are included, one also needs 

corresponding to fig. 0(b). The above integral is infrared finite because fermionic Euclidean 
frequencies po are never zero; so, unlike the bosonic case, the masses can be dropped at 
leading order mm/T . By using the same contour-trick argument that was used in Appendix 
F.l of ref. [|T|, one can easily show that 

4n = 0. (3.12) 

However, it is sometimes convenient to also know the pieces of /f^J^ corresponding to restrict- 
ing the frequency sum in various ways, and these are discussed in Appendix 

Finally, the pure gauge theory calculation required the integral corresponding to the 
bosonic piece of fig. ^(1): 

dAClg'!"^^^ ^ g -^tr [AnJ,(P)]2 , (3.13) 

where U^^, is the bosonic contribution to the vector self-energy, given by fig. [l|(a), and the 
notation 

An^,(p) = n^,(P) - n^u{o)5p^ . (3.i4) 

has been used. 

We shall need the same integral with the complete self-energy, which means we need 

dASlg'l^,^ ^ -^tr mlAP)f , (3.15) 
d^Cj^S^gH^^ ^ ijtr AnJ,(P)An;,,(P) , (3.16) 

where 11^^ is the fermionic contribution given by fig. 0(b). 

In the pure gauge theory case, all three-loop diagrams except fig. ^(1) could be reduced 
to the basketball integral /bau by the application of a few simple tricks. For instance, fig. ^(i) 
is equal to 



o o o 



(a) 



(f) 



(b) 




(g) 




(k) 



(c) 



(h) 



(1) 




(d) 



(i) 







V 



(m) 




(e) 




(j) 



FIG. 4. Diagrams contributing to the free-energy in pure gauge theory. When fermions are added, we 
include the fermionic contribution to 11^^ in diagram (1) and include the diagrams of fig. |. The crosses 
are the "thermal counter-terms" arising from the last term of ( 2.12| ), and the dashed lines are ghosts. We 
have not explicitly shown any zero-temperature counter-terms, and each diagram should be multiplied by 
the appropriate multiplicative renormalizations for vertices and propagators. 



8 



and is reduced by (1) expanding numerator factors in terms of denominator factors to cancel 
factors between numerator and denominator, such as 

P-{Q-K) = \[{K- Pf -K^-{P- Qf + . (318) 

(2) performing appropriate changes of variables to collect similar terms; and (3) using the 
identity 



(3.19) 



fp[p + QY{p + KY 2 rp{p + Q)\p + Ky 

The last identity follows by changing variables P^—P—Q—K, 

ip {P + QY{P + KY ^ ~ip {P + QY{P + KY ~ ^ ^^^ip {P + QY{P + KY ' 

(3.20) 

and then moving the first-term on the right-hand side over to the left-hand side. Unfortu- 
nately, this trick does not generalize to the case where Q+K is fermionic instead of bosonic. 
If Q+K is fermionic, then 

tp {p + qy[p + ky ^ 't{p} {p + qy[p + ky ~ ^ ^^^t{p} {P + QY{P + KY ' 

(3.21) 

and there is no simple way to solve for the bosonic integral on the left-hand side. The failure 
of this trick requires us to introduce a new fundamental integral, as was done by Parwani 
and Coriano in ref. 0:0 

" t{p}QK P^Q^K\P + QY{P + KY ■ ^^'^^^ 
If P were bosonic, this would be reducible by ( p.l9| ). 

Figs. ^ and ^ show all of the diagrams contributing to the free-energy up to three loops. 

The reductions of all the three-loop diagrams to the basic integrals are given in Appendix 0. 



^ We have adopted their notation, H^, for this integral. Their Hi and H2 correspond to our I^^^^ and I^^^ 
respectively, and their H4 is discussed in Appendix ^ 
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o e <s> Q ® 

(n) (o) (p) (q) (r) 



FIG. 5. Diagrams that must be added to fig. || to include fermions in the calculation of the free energj 

IV. RESULTS FOR INTEGRALS 



A. The fermionic basketball I^^y^ 



The derivation of I^^^i closely parallels that of the bosonic J^au in ref. [0] with the main 
difference being that the bosonic sum identity 

^g-ko|rg-|po+<?o|r- _ (cothf + |po|)e~IP"l" (4.1) 
90 

is replaced by 

^ e-ko|rg-|po+9o|r _ (cschf + |po|)e"l^°l^ (4.2) 

{90} 

where po represents bosonic frequencies and 

f = 2nTr, po = Po/27rT . (4.3) 

This has the effect of simply replacing occurrences of coth (and its small r expansion) in the 
bosonic derivation by csch (and its small r expansion). So, for instance, 

n^(^)(P) = [ d^r ^e'P-'" ( csch f--) e-l^"'!" + 0(e) , (4.4) 

where H^^^^ is the finite-temperature contribution to IT^. Similarly, 



2^4 ^oo 



327r2 Jo 



{ cschr ^ — ( — 



V rJ \ Q 

10 



(cothr — 1) 



rp4 



327r2 



/•oo / 1 \ 



C(-3) 15 5 



+ 0(6). 



(4.5) 



(47r)2 \^12^ 

The integrals in the last equation were performed using the method of Appendix |C|. Putting 
the above result together with ( |3.2| , |X^) gives 

1 



'rp2^ 2 r 



(4i)2 1,12^ 



e 47rT C(-3) C(-l) 15 5 



+ 0(e) . 
(4.6) 



The remaining terms needed to evaluate I^^^^ are discussed in Appendix 0. The final 
result for the fermionic basketball is 



'^2X2 



^'^^ (47r)2\^12; 
This agrees with the numerical result of ref. 



3 , a C'(-3) C'(-l) 173 63 



In 2 



+ 0(6). (4.7) 



B. The mixed basketball /^^n 



As has been noted by Parwani and Coriano ||^, I^l^ can be written in terms of 
J^aii by the trick of rescaling momenta {pojp) by a factor of 2: 



7-bb 
-'ball 



6{P + Q + K + R) 
PQKR P^Q^K^R^ 

2''-H i i i ^^p + Q + ^ + R) 

JP+{P} JQ+{Q} Jk+{K} Tr+{r} P^Q^K^R^ 

[-'ball ^ "-'ball ^ -'ballj ) 



(4.8) 



so that 



jbf 

-'ball 



1 /i _ 9ll-3(iA T-bb _ -'■ 7-ff 
"g l^-L ^ J -'ball g-'ball • 



To simplify notation above, we have used the short-hand 



(4.9) 



6{P + Q + K + R)^ f^-^^ 6,,+,,^ko+ro{27^y-'S'^'-'\p+q+k+r) . 



(4.10) 
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C. The new integral H3 



We now turn to the integral of ( |3.22|) , which is the one integral that is not directly 
analogous to a previous bosonic calculation. However, our attack on is inspired by our 
derivation of Jq^^ ~ |fP^^(An|^j,)^ in ref. [|I|, where we noted that the orthogonality of nj^^, 
to Pn lead to useful algebraic simplifications. We will therefore rewrite H2 in an analogous 
form. First note that is of the form 

Hs = f.^^^^,[A,{P)]\ (4.11) 

where 

Our method is to replace by something that is orthogonal to P^. So define 

/3-^^^^WP)f , (4.13) 

where 

It is easy to verify using our standard reduction tricks that P ■ J=0 and that 

H, = i/3 + + hfh -\ihi- h)lL ■ (4.15) 

Now focus on I^. The orthogonality of J to P and Lorentz invariance imply that has 
the form 

MP) = U-^P,)f{P), (4.16) 



p2 

where = (1, 0) is the four- velocity of the thermal bath. Then 



p2 y p2 
and 
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J,{P) = ^{n,-^,PA3,{P) (4.17) 



T{p} 



(4.18) 



The next simplification occurs by noting that J„ vanishes at zero temperature because 



d'^Q {2Q + P)^ 







(4.19) 



by anti-symmetry under Q —{Q+P). The large P behavior of jo{P) is therefore the large 
P behavior of its finite-temperature piece jo'^\P), which is 0(1/P^) because the 0(1/P) 
behavior of the individual pieces cancels: 



(2g + p), 



(T) 



p2 



as fermionic P — >■ oo. (4.20) 



f{p} Q'iP + Qyj 

As a result, I3 is both UV and IR finite as So we can set c? = 4 and evaluate jo for 

fermionic po in terms of the massless scalar propagator A: 

jo{P)=TY: J rfVe^P-^A(go,r-)A(po + go,rl(2go + Po) - Y^ + 0(e) 



go 
T2 



(4vr) 



T2 /• 1 ^ 



IGtt 



(cschr — cothr)e 



-|pok 



sgnpo - ^-f^ + 0{e 



8 P2 



— — / dr \ d, 



smpr 



cschr — cothr + ^rj e '^"'^sgnpo + 0{e) 



(4.21) 



57r Jo \ Pf , 

where sgnpo means the sign (±1) of po. Plugging into ( |4.18| ), doing the fermionic po sum, 
and using the identity 

1 



d^p sinpr sinps 1 
~ Air 



(27r)'^ J p^ pr ps 
to do the p integral then yields 



-9(r - s) + -9(s - r) 
r s 



(4.22) 



lz = 7; \ c/ff~^ (cschr — cothf + if ) cschr + 0(e) 

1287r"^ Jo ^ ' 

1 



9C'(-3) C'(-l) 9 9 117 ■ 



0(e). 



(4.23) 



(47r)2 \^12; 

Again, we have used the methods of Appendix |C| to do the integrals. Using ( ^.15| ) to relate 
J3 to i^a finally gives 
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(47r) 



12 



2\ 2 



3C'(-3) 3C'(-i: 



8^ + 4 ""a^ ■ 2C(-3) 2C(-1) 



9 361 57, „ 

+ -7e H \ In 2 

4 ' 160 10 



+ 0(e). 
(4.24) 



This agrees, within errors, with the numerical results of ref. 0.0 



V. RESULTS AND DISCUSSION 



The evaluation of the final basic integrals I^^^ and Jq^d closely parallel the derivation 

of 

-^qcd ^sf. and we leave the details for Appendix |F|. Combining all the results for 
individual graphs collected in Appendix our final result for the free energy is 



d^T'-< 



"9 



5 I 4rfA. 



£:) {Ca + ISf 



+ 



Air J 



Cl 



^1 38 ^-3) 148 C(-1) 64 

3 ^A-kT 3 C(-3) 3 C(-l) 5 



+ (-) C.^P 



^^n ^ l C(-3) 74C(-i: 
3 '^47rT 3 C(-3) 3 C(-l) 



1759 37, 



+ (— ^ 52F 



20, n 8C -3 16C' -1 , 1 88 , ^ 

In -^+-^ — f ^-^ — r-47E — +— ln2 

3 AtxT 3C -3 3 C -1 3 5 



An J 



105 



f24 In 2 



(5.1) 



Evaluated numerically for QCD with nf quark flavors, this is 



45 



3/2 



-|l + §nf - 0.09499 (i + A^^) + 0.12094 (l + | 
0.08662 (l + irif) In (gJl + irif) - 0.01323 (l + -|r2f) (l - |.nf) In 



+0.01733 - 0.00763 m - 0.00088 



+ 0{g' 



(5.2) 



^ We have also made a more precise numerical test of our analytic methods by computing ( 4.18 ) by brute 
force: we did the p integral and Euclidean po sum numerically and used the contour trick to get an integral 
form for jq (po ; p) i which we also evaluated numerically. The po sum converges quite quickly, and summing 
Po up to ±5/2 gave agreement with out analytic result to 0.05%. 
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FIG. 6. The dependence of the free energy F on the choice of renormahzation scale /I for 6-flavor QCD 
with as(r) = 0.1. The free energy is normahzed in units of the ideal gas result — (^oJa + -^dY-^ii^T'^ . The 
thick solid, dashed, and dotted lines are the results for F including terms through g*, 5^, and respectively. 
The light solid curve is the result plus the \ii{p,/T) term required by renormahzation group invariance. 
The light dashed curve is the result minus the g^ term. 

For QED with nf massless charged fermions with charges g^e, the free energy is 



45 



1 + |nf - 0.07916 ^ qf + 0.02328 qff^'^ 



0.00352 + 0.001341x1 ^ ) (E qjf + 0.00193 E qf 



(5.3) 
Ref. 



Our QED result agrees, within errors, with the purely numerical derivation of ref. 
also gives the QED result for the O(e^) piece. 

As in ref. |T|], we can now investigate whether the perturbative expansion of the QCD 
free energy is well-behaved for physically-realized values of couplings. Fig. |^ shows the result 
for 6-fiavor QCD when as(T)=0.1 (which corresponds to scales of order a few 100 GeV). 
The free energy is plotted vs. the choice of renormahzation scale Ji. We have taken 
1 



^ ^ ^ - /5o In ^ + §i In (1 - /5o/ (T) In ^ 



(5.4) 



where 



/3o 



(47r 



-T^ (-f Ca + |5f 



Pi 



{Att 



^ f 68^2 I 40^ c I oc 



(5.5) 
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FIG. 7. The same as fig. | but for as{T) = 0.02. 

If the expansion is well-behaved, the result for F should become more independent of 
/i as higher-order corrections are included. Instead, we see that it does not. In ref. [|l|, 
we argued that the term involves different physics than the g"^ term, and that it should 
perhaps be treated separately when discussing the behavior of the series. Ideally, one should 
calculate the free energy through g^, which is the first order that compensates for the n 
dependence of the g^ term. With our present results, however, we can at least follow ref. 
and plot results where (1) we artificially exclude the g^ term, or (2) we improve our 0{g'^) 
result with the g^ In yU term required by renormalization-group invariance: 

AF = -iAl-'y {-^ (l)" (C. + S,)hnC, - 4S.) In^j , (6.6) 

Either modification can be seen to somewhat improve the behavior of the perturbative ex- 
pansion of fig. ||. For comparison with our previous results for the pure gauge case in ref. ^ , 
figs. 1^ and ^show similar plots for the smaller couplings as{T)=0.02 and as(r)=0.001, where 
the perturbative expansion becomes progressively better behaved, as it should. For those 
readers who might be interested in the behavior of the expansion at scales of order several 
GeV, fig. ^ shows our results for the case as(T) = 0.2. 
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FIG. 8. The same as fig. | but for a^iT) = 0.001. 
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FIG. 9. The same as fig. i but for a^iT) = 0.2 and ?if = 5. 
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APPENDIX A: RESULTS FOR INDIVIDUAL GRAPHS 



Writing F=/x "^'^d^T and ignoring terms of 0(e), the diagrams of fig. ^ are given in 
Appendix A of ref. except that 

1 



in 



-{Mf + + Mi + Mi)MT , 



The diagrams of fig. 5 are given by 



(Ale) 
(All) 



c/a T{P} 



= Spg'Z'^ [{d - 2)M2h - h) + 53] , 



-^P = (-I^fOa + 52f) g\d-2) 



o -'ball ^tJball 



= S2Fg\d - 2Y + 2H, - hUi - hi 



-r = Cj,SFg\d-2)lll, 



(Aln) 
(Alo) 

(Alp) 
(Alq) 
(Air) 



The resummation of the gauge boson hne in diagram (o) is similar to the resummation of 
diagram (e), discussed in ref. In particular, 



JP{Q} 



6: 



Po 



p2 + ]Vf2 P2 



0{g\e) 



(A2) 



_Q\P + Qf Q\ 

The last equality follows from a derivation similar to the bosonic case treated in ref. [Q and 
means that /resum can be ignored at the order under consideration. 

The multiplicative renormalization constant used for the coupling is given by 

11 C^g^ , 2 S^g^ 



9ha 



Zggn' 



1 - 



+ 



+ 0{g') 



9^^ 



(A3) 



6 (47r)2e 3 (47r)2e 

The vector mass M is given by ( p.ll| ), and M| is the piece of due to the fermion 
contribution of fig. [^(b). 
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APPENDIX B: LARGE P BEHAVIOR OF U^{P) 



In Appendix B of ref. we derived the large P behavior of the bosonic U^^^^P) to be 

n^(^)(p) = 2J^^ + sjl'^ - pi) + 0{tVp') , (Bi) 



where 



The easiest way to get a similar expression for the fermionic 11^'^-^) is to note that, by scaling 
momenta (go, g) by a factor of 2, one gets 



1 _ oS-dv^ 1 



and so 



(B3) 



ii\p) = 25-°'n^(2P) - n^(P) . (B4) 

Applying this identity to ( [Bl| ) then gives the same formula, 

n^^^H^) = 2/ + 8J[^ - + OiTVP') , (B5) 

for the fermionic case but with 

Jl = (2^-^-" - l) Jl . (B6) 

APPENDIX C: NEW INTEGRALS OF HYPERBOLIC FUNCTIONS 

In Appendix C of ref. ||T|], we showed how to evaluate convergent integrals of the form 

(CI) 



Hdx Ij2 Cmnx"" coth" X + d^x^^e-^'^A 

\m,n m / 



by regulating the individual terms by introducing an extra factor of x^ and taking 5^0 at 
the end. In the current work, we need to extend our set of integrals to 
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/ = I c^npx"" coth" X csch^x + ^mX^e"'^'"^ J . (C2) 

\m,n,p m / 

To handle this, we need in addition to our previous basic regulated integrals 



oo 



dxx' = {), (C3) 

/•oo 

/ rfxx"cothx = 2-"r(2 + l)C(^ + l), (C4) 
h 

D 

dx x'e-"'' = a-^-'Til + z) , (C5) 



the new integral 



POO 

/ dx x'cschx = (2 - 2-')T(z + 1) + 1) • (C6) 
Jo 

In addition, we need to generalize our previous recursion relation to 



dx x^ coth*^ X csch^x 







oo 



dx 







and also use 



2' ft 1 

coth""^ X csch^x H -x^ coth"~^ x csch^x 



p + n — 1 p + n — 1 



(C7) 



/ (ix coth"' a; csch^x = / dx coth""*"^ a; csch^~^a; — coth*^ a; csch^~^x . (C8) 

JO JO L J 



APPENDIX D: COMPLETION OF THE CALCULATION OF I^^^-^ 

In this section, we complete the evaluation of the fermionic basketball integral. The 
derivation directly parallels the bosonic case treated in Appendix D of ref. [I]], with hy- 
perbolic cotangents becoming hyperbolic cosecants as we discussed earlier. We therefore 
refer the reader to ref. [jl| and shall here simply present the differences for a selected few 
intermediate results. 

As in the bosonic case, we write 

^^n^(T)n(o) = jf + jf + jf, (Di) 
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where 



n(o)(P) 



(47r)2e 



[nf(-)(p)-nS?)(P) 



(47r)2e Tp 



(D2) 
(D3) 



and 



^fviP) ^ 2Ji,- + (1 - 6J8JI-^^ 



p 



The zero-temperature piece 11*^'') of iP is the same as the bosonic case. One finds 



7-4 I /""^ /^pg^j^^_ 1 ^ - ''^^^ - ^ ^ Vcothf - 1) 
(47r)2 2 io f3 j 1 ^~f^6~36oJl ~2dfJ^^° ^~ 



1 /r^' 

(4vr)2 1,12, 

rp4 



+ (cschf-l + 0|+O(e) 

29 ^-3) , ^ C^(-l) 21 ^43_27 ' 
'lOC(-3) +^C(-1) 10^^+30 10 



+ 0{e) 



(4vr) 



8 4 

d- 1 



[5o(2)-rfSi(3)] 



2\ 2 



(47r)2 \, 12 
2AT 



29C'(-3) ^C'(-l) 21 43 27 



+ 5 



+ 0(e) 



4 / 47r/i2^ 



i / 

'^12^ 2 r 



(4vr)2 V12, 



(D4) 



(D5) 



(D6) 



ioc(-3) ' ^C(-i) ~ 10^'^^ 30~To^^^ 

I J^5o(l+e) + ^ [5o(2+e) - rf5i(3+e)]| 

13 /X , 2l C(-3) , 21 43 17 ■ 

20^ + To ^"i;^ + To CM) - 'cM) + ~ T + To 



(D7) 



(D8) 



where A and are defined in ref. Putting together (pq) , (|D7|) , and 



1 T 



2\ 2 



(47r)2 \,12 



1 3 , u 
H lu' ^ 



37C'(-3) ^ ^ C(-l) 301 37 



20e 10 A-nT 10 C(-3) C(-l) 120 10 

+0(e) . (D9) 



r HZ 

The result for f n(°) is given in ref. and adding it to ( f4.6| ) and ( p9D yields the final 
result (|4.7|) for the fermionic basketball. 
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APPENDIX E: THE PIECES OF 1*^^ 



As mentioned earlier, the fermionic setting-sun integral I^^^ vanishes when particle 
masses are ignored. However, it is useful to also know the piece corresponding to 



This can be easily evaluated by relating it to the comparable bosonic piece using 
then scaling momenta by 2: 



-^n^(p) = [2'-''u^{2P) - n^(P) 

6„ 



S„ (1 - 6. 



(2- - 1) 



T2 



PQ P^Q\P + Qf 
ln2 + 0(e). 



where the last equality follows from the bosonic result of ref. [||] that 



4^ 



T2 



fP^Q^{P + Q)^ (47r)2 



1 2T 1 

h In 

4e ft 2 



+ 0(m,e) 



(El) 



and 



(E2) 



(E3) 



APPENDIX F: DERIVATIONS OF ij^ AND 



The calculations of J^^j and I^^^ also closely parallel those of the purely bosonic case 
in Appendix H of ref. |^. In that previous work, we found it algebraically convenient to 
rewrite the bosonic contribution IlJ^^ of fig. |l|(a) to the vector self-energy in terms of 



n 



(Fl) 



via 



d-2^ 



n^,(P)-2(PV-P^P, 



QQ^{p + Qy 



(F2) 



23 



In order to make our presentation of the fermionic case as similar as possible to the bosonic 
one, we shall do the same for the fermionic contribution 11^^ of fig. 0(b) :0 

1 



{Q} Q\p + Qy 



(F3) 



where 



TTf (2g + P),{2Q + P\ 



(F4) 



And now let's focus on computing 



sqed 



mup)f 



sqed 



AW(P)AW{P) 



(F5) 



IP P4 ' sqea ^ p4 

Following through the same steps as in Appendix H of ref. ||^, one gets the obvious gener- 
alizations of eqs. (H27-29) of that reference: 



t K^J)(P)]^ = [nf(^)(P)]% Aid - 2)h,!l + 0(e) , 

;^Anjr(^)An;(J)(P) = nb(^)(P)nf(^)(P) + Aid - 2)b,bJ, + 0(e) , 

i ^^ll\P)mP) = n^^'H^)n^°HP) + 2|^/i^^ ^n(o)(p) , 

t ^nJFHP)ni°j(P) = nM-)(P)n(o)(P) + i,n(o)(P) 



P P4 



d 



l^^Jn(o)(P) 



(F6) 
(F7) 

(F8) 
(F9) 

(FIO) 



Summing these results, and incorporating the results for the assorted basic integrals, gives 



1 (T^^^' 



11 , u IC -3 20 C -1 281 

— + 11 n-^ + — f + — Z + 4^ + 13 n2 

6e 47rT 3 C -3 3 C -1 60 



rbf 

sqed 



(47r)2 \^12^ 



+0(e), 

59 _ ^ + 11 C'(-3) 70 C'(-l) 



(Fll) 



12e 2 A-kT 6 C(-3) 3 C(-l 
+0(e) . 



2063 169, „ 
?7e - — r- + — In 2 



120 10 



(F12) 



^ The main advantage to this is simphcity of presentation; one could just as easily do the calculation with 
nf directly. 
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Using ( |F2|) and ( p3| ) and our standard reduction tricks yields:^ 



4d = <ed + 4(rf - 3)/« u - 16{d - 2)h^ -ln^(P) 



1 [T^V 



(47r)2 \^12^ 



40 , u 32 C' -3 224C' -1 124 448, 

— + 80 In — ^ — f + ^ — f + 167e + + In 2 

3e A-kT 3 C -3 3 C -1 '35 



+0(e) , (F13) 

= (2 - rf)/.^'ed - 3(rf - 2)d + 8(rf - 2)hf^ ^^\P) + 2(^; - 2fb{^^ ^^\P) 

Ji 38C'(-3) 104C'(-1) 



(47r)2 I 12 J 



29 

'37 ''47rT 3 C(-3) ' 3 C(-l) 



251 398 
■ 167e - — + — In 2 - 96 ln(27r) 
10 5 



+ 0(e) . 



(F14) 



Before leaving this section, we note that one can use our results to derive the basic 
integral 

{Q ■ Kf 



Ha = 



(F15) 



fp{QK} P^Q-'K^P + Q)\P + KY 

defined by Parwani and Coriano. By applying our usual reduction techniques to If^^edi 
using the fact that Jf^^ = 0, one can derive that 



(F16) 



Solving for and using our results for J^^j and J^au then yields 



2 r 



(47r)2 \^12 



5 5 J 

2At^A^W 6 C(-3) ' 6 C(-l 



lC'(-3) , 7C'(-1) ^1 ^23 
+ r^+24 



-In 2 

5 



0(e), 
(FIT) 



which agrees with the numerical result of Parwani and Coriano within errors. 



* Because of our slightly different methods of bookkeeping of infrared divergences, our reduction of I^^^ 
differs from a similar reduction in ref. |^ by Ml^{l — (5p^)P^^n^(P). There is a canceling difference in our 
treatment of diagram (o) of fig. ||. 



25 



REFERENCES 



[1] P. Arnold and C. Zhai, University of Washington preprint UW/PT-94-03, [hep- 
I ph/94082T6 . 



[2] C. Coriano and R. Parwani, Argonne National Lab preprint ANL-HEP-PR-94-02 (re- 
vised) (1994), |hep-ph/9405343| ; R. Parwani, Saclay preprint SPhT/94-065 (1994), [h^ 
ph/9406318| ; R. Parwani and C. Coriano, Argonne National Lab preprint ANL-HEP-PR- 
94-32 (1994), |h^li/940926£ . 



26 



